Two black hole initial data 
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The Misner initial data are a standard example of time-symmetric initial data with two apparent 
horizons. Compact formulae describing such data are presented in the cases of equal or nonequal 
masses (i.e. isometric or non-isometric horizons). The interaction energy in the "Schwarzschild + 
test particle" limit of the Misner data is analyzed. 

PACS numbers: 04.20.Ex, 04.70.Bw 
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I. INTRODUCTION 

In p| Misner proposed two-body wormhole initial data and in 0] he used the method of images to describe the 
time-symmetric initial data for an arbitrary number N of particles (see also for a review of time-symmetric initial 
data). Such data may be viewed as a collection of Einstein-Rosen bridges connecting two isometric sheets. If we 
restrict ourselves to just one of them, we can view such initial data as a three-dimensional metric on a manifold 
with internal boundary. The boundary consists of N minimal surfaces representing surfaces of the "black holes" 
corresponding to the particles. In this paper I am concerned with the N = 2 case only. If the masses of black holes 
are equal (i.e. when the horizons are isometric) then there is a well-known formula for the three-metric: 



$ 4 ( d/j 2 + dif + sin 2 r] dip 2 ) , 



where 



Vd 



n£Z y/cosh(fi + 2nfi Q ) - cos 77 
Here (p, 77, tp) are the bispherical coordinates defined by 
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cosh /1 — cos 7] 
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cosh /i — cos 77 ' 
sinh /i 



cosh /1 — cos 77 



(3) 



There are two parameters in (j2J: d and jiQ. The jj, coordinate ranges from — /j,q to /_tn- The /1 = ±/xq surfaces are 
minimal, which translates to the Neumann boundary condition for $: ^$^=±^0 = 0. The conformal factor <f> 
satisfies an elliptic equation 



sin 77 ( A - - I <I> 



(4) 



where A is the Beltrami-Laplace operator associated with the metric d/j 2 + dr/ 2 + sin 2 77 dip 2 and sin ij is the volume 
element. 

It is obvious that the metric g is conformally related to the Euclidean metric on R 3 minus two balls which correspond 
to > fiQ. This transformation reads: 



,g = * 4 d 2 (d^+ dy z + dz z ) , 



(5) 
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where the conformal factor \I/ is related to $ by the formula ^> — £f _1 / 2 $-y/cosh/x — cos 77. The yd factor is extracted 
from $ in order to normalize the value of W at infinity to 1. Now the parameter 2d has the interpretation of the 
Euclidean distance (with respect to the metric d 2 (dx 2 + dy 2 + dz 2 )) between the points corresponding to /i = ±00 
(* = ±1). 

To generalize the metric g to the case of nonequal masses one takes fi G [—a, b] with a, b > 0, imposes the Neumann 
boundary conditions and solves equation Q for <£>. Such metrics are well known and used to construct initial data. 
There are effective formulae which give the conformal factor, see for example |4j. However, it seems that no closedform 
formula has been published, for example the formulae from 2] include series of operators of inversion with respect to 
spheres in Euclidean space and the formulae in £| are recursive. Below I present a compact form of the metric. First I 
derive a new version of formula (J2J) in the equal-mass case, see (JZJ and JSJ. Then I present the corresponding formulae 
I|1L)| I and (|lll) in the case of nonequal masses. I also use (|11|) to calculate the interaction energy in the test-body limit 
of the Misner initial data, that means when one mass is much smaller than the other. 



II. MAIN RESULT 



Let us begin with 

1 00 1 



1 1 1=0 > 



where r< := min(|r|, |r'|), r> := max(|r|, |r'|), (3 — the angle between r, r'. Taking f = e z and \r'\ = exp(a), we get 



\e z — r'\ — V2e a / 2 \J cosh a — cos 77 , 

where rj is a spherical angle of the vector r'. Hence we have the following formula for terms in (J2J: 

1 00 

v /2^ e -l«K'+ 1 /2)p i ( cosr? ) ) (6 ) 



v^osh a — cos r\ 



1=0 



where a ^ 0. If we expand the summands in (J2J) using (j^Jl and sum the n-indexed series then we arrive at the conformal 
factor for the equal-mass case rewritten as: 

^ sinh[/j (Z + 1/2)] 



the formula being valid for // ^ 0. We can also extract the singular term ^ cosi f^_ cosri an d get the following version of 
0: 



*= , f (8) 
Vcosh [J, — cos T] 



K ' sinhM* + l/2)] ' 



with the series being uniformly convergent for jj, £ [—Ho, A*o]- 

In the case of non-equal masses the conformal factor is given by the following formula: 

\fd °° 

— 2na—2(n—l)b T" ty — 2na— 2nb (n— l)a-\-2nb + Q2 na-\-2nb ) , (9) 

V cosh /i — cos 77 ^ 

where 

1 

Q a 



^cosh a — cos rj 



Simple symmetry argument shows that such $ indeed satisfies the Neumann boundary conditions at /i = —a, /i = b. 
As before we use (jSJl and perform summation over n to get $ in the following compact form: 

d> = ^vp i( co S ,) cosh ^ - b + w + + b - \m ± v 2 )] , do) 

^ smh[(a + b)(l+ 1/2)] 
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or, if we extract the singular term, 

yd 

*= , = =+ (11) 
y cosh — cos 77 

~ e -M(' + l/2)( e 26(/+l/2) + 1) + e M(i+l/2) ( e 2o(J+l/2) + -y\ 

V2i2^fl(cOS»j) e2 (a +b )( i + l/2) _ 1 ■ 

Again formula l|10|) is valid for /i ^ and the series in converges uniformly for \i e [—a, 6]. If we substitute 
a = b = fio then ifTUI) and (fTTft reduce to and ©, respectively. 

The main result here are formulae (|10fl and It is worth noting that the conformal factor is expanded in a 

series of orthogonal polynomials, which is a desirable feature for numerical treatment. 

III. SCHWARZSCHILD + TEST-BODY LIMIT 

Let us now analyze the "Schwarzschild + test body" limit of Misner initial data. First, we observe that if we pass 
to the limit b — > 00 in l|ll(l then we get Schwarzschild initial data in bispherical coordinates. The /i = —a sphere 
represents the minimal surface in Schwarzschild initial data and /j, = 00 becomes a regular point. The mass of such 
data is 

2d , s 

(12) 



sinh a 

Second, we treat e := exp(— b) as a small parameter, that means s <C exp(— a), and expand (|11|) in e. We get the 
following formula for the conformal factor: 

$ = / f < 13 > 
Vcosh/i — cos 77 

°^ e - M (*+l/2^ 



Z-^i e 2a(l+l/ 
1=0 

Let us denote by mi and m-i the masses of the (i = —a and the /1 = b surfaces respectively. By "mass of a surface" I 
mean here the square root of its area divided by \/16ir. Using (|13fl we find the individual masses of the surfaces and 
the ADM mass M of the whole system: 

mi = + 8d{e- a + e- 2a )e + 0(e 2 ) , (14) 

sinh a 

m 2 = Ad(l + e - a ) 2 e + 0{e 2 ) , (15) 

M = mi + to 2 - 8d(e~ a + e - 2a )e + 0(e 2 ) . (16) 

In the Newtonian gravity the interaction energy equals E = — JT£ ™g e and for the Misner data the interaction energy 
should converge to this value in the Newtonian limit, that means when the distance increases to infinity. We would 
like to write the interaction energy in the form E = — J™l^. c l, where 7 is a function of distance and converges to 
1 as the distance goes to infinity. The two most obvious choices for the distance parameter are D := 2d and the 
length of the 77 = n geodesic connecting two horizons which we denote by L. The distance D is quite close to the 
distance studied by Brill and Lindquist 0, in the sense that it is the Euclidean distance between two points lying 
inside minimal surfaces. In this case we get the following formula for the interaction energy: 

M-mi- m 2 = fl ^ ) + 0(e 2 ) . (17) 

D V i+v 1 + (^) I 

Passing now to the second case, we can calculate the geodesic distance L as the integral 

1 



L = d I {y(x = 0,y = 0,z)) 2 dz + O(e), 

Zl 
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where z\ 



' cxSfoTT ^ s * ne z coor dinate 10 of the point (/x = —a, n = ir). The result is 

1 1 



2d 



III T 



1 _ £i nil 1q 1 - zp 

2 2 2<i zi — zq 8d 2 \z\ — z$ 1 — zq 



0(e) 



with zq = — ° being the z coordinate of the center of the removed ball corresponding to the /i 



surface. As the result we get the interaction energy equal to 

toito 2 

Ai — mi — TO 2 — 



7l + 0(e 2 ) 



minimal 



(18) 



The factor 



1l 



L 

2d. 



II 1 1 

2d 



i + V 1 + (W 



can be expressed in zeroth order in e (using 1)14(1) as a rather complicated function of a single parameter — . In the 
d — > 00 limit (L — > 00) it converges to 1, giving the correct Newtonian limit of (|18|) . 
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